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METHOD AND SYSTEM TO SCALE DOWN A DECISION TREE-BASED 
HIDDEN MARKOV MODEL (HMM) FOR SPEECH RECOGNITION 

FIELD OF THE INVENTION 

The present invention relates generally to speech processing and to automatic 
speech recognition systems (ASR). More particularly, the present invention relates to a 
method and system to scale dovm a decision tree-based hidden markov model (HNDvI) 
for speech recogrution. 

BACKGROUND OF THE INVENTION 

A speech recognition system recogrdzes a collection of spoken words ("speech") 
into recognized phrases or sentences. A spoken word typically includes one or more 
phones or phonemes, which are distinct sounds of a spoken word. Thus, to recognize 
speech, a speech recogrution system must determine relationships between the words 
in the speech. A common way of determining relationships between words in 
recognizing speech is by using a general-purpose acoustical model ("general model") 
based on a hidden markov model (HMM). 

Typically, a HMM is a decision tree-based model in which the HMM uses a 
series of transitions from state to state to model a letter, a word, or a sentence. Each 
arc of the transitions has an associated probability, which gives the probability of the 
transition from one state to the next at an end of an observation frame. As such, an 
unknown speech signal can be represented by ordered states with a given probability. 
Moreover, words in an unknown speech signal can be recognized by using the ordered 
states of the HMM. The HMM, however, can place a heavy burden on system 
resources. 

Thus, a challenge for speech recognition systems is how to utilize system 
resources for improving the performance of using a general model such as the HIvI^/I. 
A disadvantage of using the general model is that it is trained for broad use from a 
very large vocabulary, which can lead to poor performance for special applicatior\s 
related to specific vocabulary. For example, a mismatch between speaker 
characteristics, transmission chaimels, training data, etc., can degrade speech 
recognition performance using the general model. Another disadvantage of the 
general model is that it requires extensive computation costs and high resource 
utilization. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The features and advantages of the present invention are illustrated by way of 
example and not intended to be limited by the figures of the accompanying drawings 
in which like references indicate similar elements, and in which: 

FIG. 1 is a diagram illustrating an exemplary digital processing system for 
implementing the present invention; 

FIG. 2 is a flow chart illustrating an operation to scale a general model into a task 
specific model according to one embodiment; 

FIGS. 3 A through 3D are exemplary diagrams illustrating how decision trees can 
be scaled; 

FIG. 4 is a flow chart illustrating a trim-down operation according to one 
embodiment; 

FIG. 5 is a flow chart illustrating a task adapting and interpolating operation 
according to one embodiment; and 

FIG. 6 is a flow diagram illustrating an interpolation process according to one 
embodiment. 

DETAILED DESCRIPTION 

According to one aspect of the invention, a method and system are provided in 
which a decision tree-based model ("general model") is scaled down ("trim-down") for 
a given task. The trim-down model can be adapted for the given task using task 
specific data. The general model can be based on a hidden markov model (HMM) and 
can be trained through a large scale general purpose (task-independent) speech 
database. 

By allowing a decision tree-based acoustic model ("general model") to be scaled 
according to the vocabulary of the given task, the general model can be configured 
dynamically into a trim-down model, which can be used to improve speech 
recognition performance and reduce system resource utilization. Furthermore, the 
trim-down model can be adapted/adjusted according to task specific data, e.g., task 
vocabulary, model size, or other like task specific data. 

FIG. 1 is a diagram illustrating an exemplary digital processing system 100 for 
implementing the present invention. The speech processing and speech recognition 
techniques described herein can be implemented and utilized within digital processing 
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system 100, which can represent a general purpose computer, portable computer, 
hand-held electronic device, or other like device. The components of digital 
processing system 100 are exemplary in which one or more components can be 
omitted or added. For example, one or more memory devices can be utilized for 
digital processing system 100. 

Referring to FIG. 1, digital processing system 100 includes a central processing 
unit 102 and a signal processor 103 coupled to a display circuit 105, main memory 104, 
static memory 106, and mass storage device 107 via bus 101. Digital processing system 
100 can also be coupled to a display 121, keypad input 122, cursor control 123, hard 
copy device 124, input/output (I/O) devices 125, and audio/speech device 126 via 
bus 101. 

Bus 101 is a standard system bus for communicating information and signals. 
CPU 102 and signal processor 103 are processing units for digital processing system 
100. CPU 102 or signal processor 103 or both can be used to process information . 
and/or signals for digital processing system 100. Signal processor 103 can be used to 
process speech or audio information and signals for speech processing and 
recognition. Alternatively, CPU 102 can be used to process speech or audio 
information and signals for speech processing or recognition. CPU 102 includes a 
control unit 131, an arithmetic logic unit (ALU) 132, and several registers 133, which 
are used to process information and signals. Signal processor 103 can also include 
similar components as CPU 102. 

Main memory 104 can be, e.g., a random access memory (RAM) or some other 
dynamic storage device, for storing information or instructions (program code), which 
are used by CPU 102 or signal processor 103. For example, main memory 104 may 
store speech or audio information and iristructions to be executed by signal processor 
103 to process the speech or audio information. Main memory 104 may also store 
temporary variables or other intermediate information during execution of 
instructions by CPU 102 or signal processor 103. Static memory 106, can be, e.g., a 
• read only memory (ROM) and/or other static storage devices, for storing information 
or instructions, which can also be used by CPU 102 or signal processor 103. Mass 
storage device 107 can be, e.g., a hard or floppy disk drive or optical disk drive, for 
storing information or instructions for digital processing system 100. 



Display 121 can be, e.g., a cathode ray tube (CRT) or liquid crystal display (LCD). 
Display device 121 displays information or graphics to a user. Digital processing 
system 101 can interface with display 121 via display circuit 105. Keypad input 122 is a 
alphanumeric input device for communicating information and command selections to 
digital processing system 100. Cursor control 123 can be, e.g., a mouse, a trackball, or 
cursor direction keys, for controlling movement of an object on display 121. Hard 
copy device 124 can be, e.g., a laser printer, for printing information on paper, film, or 
some other like medium. A number of input/output devices 125 can be coupled to 
digital processing system 100. For example, a speaker can be coupled to digital 
processing system 100. Audio/speech device 126 can be, e.g., a microphone with an 
analog to digital converter/for capturing sounds of speech in an analog form and 
transforming the sounds into digital form, which can be used by signal processor 203 
and/or CPU 102, for speech processing or recognition. . 

The speech processing techniques described herein can be implemented by , 
hardware and/ or software contained within digital processing system 100. For 
example, CPU 102 or signal processor can execute code or instructions stored in a 
machine-readable medium, e.g., main memory 104, to process or to recognize speech. 

The machine-readable mediiom may include a mechanism that provides (i.e., 
stores and /or transmits) information in a form readable by a machine such as 
computer or digital processing device. For example, a machine-readable medium may 
include a read only memory (ROM), random access memory (RAM), magnetic disk 
storage media, optical storage media, flash memory devices. The code or instructions 
can be represented by carrier wave signals, infrared signals, digital signals, and by 
other like signals. 

FIG. 2 is a flow chart illustrating an operation 200 to scale a general model into a 
task specific model according to one embodiment. Operation 200 includes two parts. 
The first part relates to how the general model is scaled by the vocabulary of the given 
task to create a trim-down model. The second part relates to adapting the trim-down 
model based on task specific data to create a task specific model. These two parts of 
operation 200 can be used separately or in combination. 

Referring to FIG. 2, at operations 202 and 204, vocabulary knowledge of a given 
task and a creneral model are input to a trim-down operation 206 for reducing the 
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general model according to the vocabulary of a given or specific task into a "trim- 
down model". 

The general model may be a general-purpose acoustical model such as the HMJvI. 
For example, the general model may be a decision tree-based HIvI>/L The general 
model may include an exhaustive list of words for a given language. The vocabulary 
knowledge of the given task, however, may include only a small subset of the general 
model. For example, the vocabulary knowledge of the given task may be related to 
"weather reporting" and include words and phrases such as, for example, "cloudy," 
"rain today," "temperature," "hunudity," "the temperature today will be," etc. 

At operation .206, the general model is trimmed down according to the 
vocabulary of the given task to create a "trim-down model." The trim-down operation 
206 allows a general model to be scaled for use by specific vocabulary instead of 
having to use the exhaustive vocabulary of the general model. 

At operation 208, the trim-down model is adapted with task specific data to 
configure the trim-down model according to specific task requirements such as, for 
example, task specific vocabulary, model size, etc. For example, the trim-down model 
can be trained for a task dependent model- 

At operation 210, a task specific model can be obtained after adapting the trim- 
down model. For example, a general model can be scaled down for a given task such 
as, for example, weather reporting and adapted with weather reporting specific data to 
create a task specific model for weather reporting. Thus, operation 200 provides a fast 
and accurate task specific model for recognizing speech of the given task by 
downsizing or scaling a general model for the given task. 

FIGS. 3A through 3D are exemplary diagrams illustrating how decision trees 
can be scaled. For purposes of illustration, the decision trees are trees based on 
decision tree state clustering for a general model such as the HMM. In the following 
examples of FIGS. 3A through 3D, child nodes or leaf nodes of the same parent node 
represent HMM states of similar acoustic realization. Parameters under child or leaf 
nodes are clustered down to a smaller size. Each HMM state represents a triphone 
unit. A triphone unit includes a phone or phoneme having a left and right context. 

Furthermore, each node of the decision trees is related to a question regarding 
the phonetic context of the triphone clusters with the same central phone. A set of 
states can be recursively partitioned into subsets according to the phonetic questions at 
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each tree node if traversing the decision tree from the root to its leaves. States reaching 
the same child or leaf node on the decision tree are regarded as similar and the number . 
of child or leaf nodes determines the number of states. 

The following exemplary decision trees illustrate the rxiles for merging nodes of 
a decision tree of a general model. That is, if all the child or leaf nodes from a parent 
node do not occur in the given task, the child or leaf nodes will merge with the parent 
node. In addition, if all the child or leaf nodes from a -parent node do occur in the 
given task, the child or leaf nodes will not merge with the parent node. Furthermore, 
if any parent node having incomplete descendants (i.e., descendants that do not occur 
in the given task), it will be replaced by its child or leaf node, which does have 
complete descendants. 

FIG. 3A and 3B illustrate if all of the child or leaf nodes from a parent node do 
not occur in the given task, the child or leaf nodes will merge with the parent node. 
Referring to FIG. 3A, an exemplary decision tree is shown having a plurality of nodes 
302 through 310. A node labeled with a line (' through it refers to a node that 
occurs in a given task. Node 308 refers to a node that occurs in the given task. For 
example, node 308 may refer to the word "cloudy" for the given task of "weather 
reporting." For one implementation, during a trim-down operation, nodes like node 
308 and 310 will merge with its parent node 304 because not all of the child or leaf 
nodes from its parent node 304 occur in the given task. Referring to FIG. 3B, another 
exemplary decision tree is shown having a plurality of nodes 312 through 320. The 
decision tree has no nodes labeled with a line ("-") through it. Thus, for another 
implementation, during a trim-down operation, child or leaf nodes like nodes 318 and 
320 will merge with its parent node 314 because they do not occur in the given task. 

FIG. 3C illustrates if all of the child or leaf nodes from a parent node do occur in 
the given task, the child or leaf nodes will not merge with the parent node. Referring 
to FIG. 3C, an exemplary decision tree is shown having a plurality of nodes 322 
through 330. A node labeled with a line ("-") through it refers to a node that occurs in 
a given task. The nodes 328 and 330 refer to a node that occurs in the given task. 
Thus, for another implementation, during a trim-down operation, nodes 328 and 330 
will not merge with its parent node 324. 

FIG. 3D illustrates that if any parent node having incomplete descendants (i.e., 
descendants that do not ocoir in the given task), it will be replaced by its child or leaf 
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node, which does have complete descendants. Referring to FIG. 3D, an exemplary 
decision tree is shown having a plurality of nodes 332 through 342. The parent node 
333 has child or leaf nodes that are incomplete, however, child or leaf node 338 has 
complete descendants, i.e., nodes 340 and 342. Thus, node 338 will merge with parent 
node 333. Node 336 may also merge with parent node 333 following the rules set forth 
in FIGS. 3A and 3B. 

FIG. 4 is a flow chart illustrating in further detail the trim-down operation 206 
of FIG. 2 according to one embodiment. Referring to FIG. 4, at operation 402, a 
decision tree of a general model is obtained. The decision tree can be similar to the 
exemplary decision trees shown in FIGS. 3A through 3D. 

At operation 404, leaf nodes of the are labeled if they occur in the given task 
vocabulary. For example, as shown in FIG. 3A, node 308 is labeled as having occurred 
in the given task. 

At operation 406, leaf nodes are merged if necessary. For one implementation, 
as illustrated in FIGS. 3 A and 3B, if all of the child or leaf nodes from a parent node do 
not occur in the given task, the child or leaf nodes will merge with the parent node. 
For another implementation, as illustrated in FIG. 3C, if all of the child or leaf nodes 
from a parent node do occur in the given task, the child or leaf nodes will not merge 
with the parent node. For another implementation, as illustrated in FIG. 3D, if any 
parent node having incomplete descendants (i.e., descendants that do not occur in the 
given task, it will be replaced by its child or leaf node, which does have complete 
descendants. 

The merged nodes in the trim-down operation 206 can be represented by mixed 
gaussian distributions to achieve high performance in speech recognition systems. For 
example, given two leaf ("child") nodes n^ and n^, with Gaussian distributions 
Gj = N(/x, , ) and = N{]U.^ respectively, using a Baum-Welch ML estimate, 
the following two equations can be used: 



where X = {speech data aligned to Gaussian G, with occupancy count /(.x) for 
each data x}, a = ^Yix) is total occupancy of Gaussian G, in the training data. 



Furthermore, assuming that both sets of data X and Y are modeled by the combined 
Gaussian G = l^iM^^)' i-^-/ when the two leaf nodes n, and ri, are merged together. 
The following merged mean and variance can be computed as follows: 
a b 

= 7 1^1 + 7^2 W 

where b = ^/(y) is total occupancy of Gaussian G^ in the training data. 

FIG. 5 is a flow chart illustrating a task adapting and interpolating operation 500 
according to one embodiment. Referring to FIG. 5, at operation 502, a trim-down 
model of a general model is obtained- For example, the operation 206 can be used to 
obtain a trim-down model. 

At operation 504, the trim-down model is trained with task specific adaptation 
data. The task adaptation data can be used to adapt the general model (using the trim- 
down model) for the vocabulary o the given task. At operation 506, a task dependent 
model is derived after the training process of operation 504 using the trim-down 
model and the task specific adaptation data. 

At operation 508, the trim-down model and task dependent model can be 
interpolated. The interpolation process can use an approximate maximum a posterior 
(AMAP) estimation process as will be explained in FIG. 6. 

At operation 510, a specific task model is generated after interpolating the trim- 
down model and the task dependent model. 

FIG. 6 is a flow diagram illustrating an interpolation process (AMAP) 600 
according to one embodiment. Referring to FIG. 6, at functional block 602, adaptation 
data is input to a Baum-Welch process to obtain a posterior probability. The posterior 
probability can be obtain as follows. 

Suppose that y,(5,) is the probability of being at state s, at time t given the 
current HMM parameters, and 0,,(j,) is the posterior probability of the ith mixture 
component. 
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(I>uis,) = Pico, |X,.5,) 

Xpi(Oj\x,)N{x,;fij,fij) 

At functional blocks 606, 608, and 610, an approximate N^IAP (AMAP) 
estimation scheme can be implemented by linearly combining the general model and a 
task dependent model (TD model) with task specific counts for each component 
density. 

=A(.),."-.(1-A)(.x> - 
n,^ =A«/'+(l-A)n/° 

Where the superscripts on the right-hand side denote the data over which the 
following statistics (counts) are collected during one iteration of the forward-backward 
algorithm. 

{x).^^^r^is,)<t>M)^x, 



Where the weight A controls the adaptation rate. Using the combined counts, 
we can compute the AMAP estimates of the means and covariances of each Gaussian 
component density from 



/ N AMAP 
AMAP _ \-^/, 
Ml — AMAP 



T\ AMAP 



AMAP \XX )i AMAPr AMAP^T 

= JmP ^f^i ^ 

At functional block 612, a task adaptation model or a task specific model is 
obtained by the combination of the general model and the task dependent model. 

The following is exemplary high level code for implementing a trim-down 
operation to scale a general model. For example, the trim-down operation as 
described herein can be implemented using the following exemplary code. 
Furthermore, the trim-down operation can be implemented with any standard 
programming language such ais, for example, "C" or "C++." 
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^include "hmm/hmmphys.h" 
delude "hmm/hmmparam.h" 
rfinclude "hmm/hnmninap.h'* 
rtadude "himn/hmmlist.h" 
^include "hmm/phnench" 
#include "hmm/ phnname.h" 
rrinciude "hmm/mono.h" 
^include "dcstree/question-h" 
#include "statis/statis.h" 
rtoclude "dcstree.h" 
rrinciude "dsizemodel_voc.h" 
delude "logmathJi" 
^include "exscan.h" 

extern HMMPhysical *hminphys; 
extern HM^/tParam *hinmparain; 
extern HMIvIMap *hminmap; 
extern DcsTreeQuestion *dtq; 
extern STATIS *hinmocc; 
extern ExStateScan *exscan; 
extern MonoPhone *mono; 
extern float ^newWgtSet; 
extern float *newMeanSet; 
extern float *newVarSet; 
extern int *newNmixSet; 
extern int *newMixSet; 

/* This function combine the nodes with none or only one child node can be traversed from the 
vocabulary list */ 

int DSizeModel::VocMergeNodes^M(aDcsTreeNode ♦node^int magic,int newstate) { 
printfC'entering node%d\n",node->id); 
node->magic=magic; 
if (node->Voc_Exist>=0) return newstate; 
if (node->Voc_Exist==-l) { 
if (!node->yes) return newstate; 
if (node->yes->Voc_Exist==-l) 
newstate=VocMergeNodes_M(node->yes,magic^ewstate); • 

if (node->no->Voc_Exist==-l) 
newstate- VocMergeNodes_M(node->no,magic,newstate); 

if ((node->yes->Voc^Exist>=2) | | (node->no->Voc_Exist>=2)) { 
node->Voc,Exist=node->yes->Voc.Exist+node->no->Voc_Exist; 

return newstate; 

if ((node->yes->Voc_Exist==l)&&(node->no->Voc_Exist==l)) { 
node->Voc_Exist=2; 
return newstate; 

if ((node->yes->VcK:_Exist>=0)«k&(node->no->Voc_Exist>=0)) { 
float cost=MergeGaiassianCost(node->yes,node->no); 
//if (cost>=0.8e+10) { printfrcost=%13.2f\n");retum newstate;} 

printfC'Merge node%d and node%d, the lost cost=%10.2f\n",node->yes->id,node.>no->id,cost); 

for (int sss=0; sss<anux; sss++) { 
FinalCGLTopl2[sssl.combCost=CGLTopl2[sss].combCost; 
FinalCGLTopl2[sss].combWeight=CGLTopl2[sss].combWeight; 

for (int ttt=0;ttt<DSvecSize;ttt++) { 
FinalCGLTopl2[sss].combMean[ttt]=CGLTopl2[sss].combMean[ttt]; 
FinalCGLTopl2[s5sl.combVar(tttl=CGLTopl2[sssl.combVar[tttl; 

} 

) 
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) 

CopyParaml(node->yes,node->no,newstate); 

printffnew state %d is generated successhilly!\n",newstate); 

node->Voc_Exist=node->yes->Voc_Exist+node->no->Voc,Exist; 

node->state=newstate+newstatebase; 

node->occ=aode->yes->occ+node->no->occ; 

newstate-M-; 

return newstate; 

I 

/* This function mark the nodes with one child node's Voc_Exist=0 */ 

aDcsTreeNode * DSizeModel::VocMergeNodes.Refine(aDcsTreeNode *node,int magic,aDcsTreeNode 
*teinp_node) { 

printfC'entering node%d\n",node->id); 

node->mag;ic=magic; 

if (node->Voc_Exist<=l) return temp_node; 
if (node->Voc_Exist>=2) { 

if (!node->yes) { printf("Error\n"); return temp.node;} 

temp_node=VcK:MergeNodes_Refine(node->yes,magic,temp_node); 

temp_node=VocMergeNodes_Refine(node->no,magic,temp_node); 

if (((node->yes->Voc^Exist>=2)&&(node->nch>Voc_Exist==0)) | | ((node->no- 
>Voc_Exist>=2)&&(node->yes->Voc_Exist==0))) { 
if (temp_node) { 

printf("Starting refinement,%d node is being replaced by %d node\n", node->id, temp_node->id) 
node->occ=node->yes->occ+node->no->occ; 

node->Voc_Exist=1000; /* mark with a large number which can not be achieved in general cases ' 

return node; 

.. } ■ 

return temp^node; 

if ((node->yes->Voc_Exist==l)&&(node->no->Voc.Exist==l)) { 
return node; 

} 

return node; 

} 

float DSizeModel::MergeGaussianCost(aDcsTreeNode *a,aDcsTreeNode *b) { 
int esidxl=a->state; 
int esidx2=b->state; 

float al,bl,ininninix,maxranix; 
irit mixl,mix2; 
int nmixl,nmix2; 

float *meanl,*varl/mean2/var2,*meanl.r,*varl_r/mean2_r,*var2.r; 
int *newnmix=newNmixSet; 
int *newmix=newMixSet; 
if (esidxl<newstatebase) { 

if (a->occ<=0.0) a->occ=hmmocc->GetOcct(esidxl); 

mixl=hmmparam->ExState(esidxl)->mix; 

nmixl=hmmparam->ExState(esidxlV>nmLx; 

meanl=hmmparam->Mean(mixl); 

varl=hmmparam->Var(mixl); 
}else{ 

nrrdxl=newnmix[esidxl-newstatebasel; 

meanl.r=newMeanS€t+newmix[esidxl-newstatebase]*DSvecSize; 
varl_r=newVarS€t+newmix[esidxl-newstatebasel*(DSvecSize-5-l); 

} 

if (esidx2<newstatebase) { 
if (b->occ<=0.0) b->occ=hmmocc->GetOcct(esidx2); 
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rTux2=hmmparam->ExState(esidx2)->mix; 
ruTux2=hmmparam->ExState(esidx2)->riinix; 
mean2=hminparam-> Mean( tnix2); 
var2=hnunparam->Var(nux2); 
lelsel 

njrux2=newnirdx[esidx2-newstatebasel; 

mean2_r=newMe2mSet+newmix[esidx2-newstateba5el*DSvecSize; 
var2_r=newVarSet+newmix[esidx2-newstatebase]*(DSvecSize-rl); 

I 

float *wgtl,*wgt2; 
if (esidxl<newstatebase) 
wgtl=hmmparam->Weight(hininparam->ExState(esidxl)); 

else 

wgtl=newWgtS€t+newtnix(esidxl-newstatebase]; 

if (esidx2<newstatebase) 
wgt2=hininparam->Weight{hinmparam->ExState(esidx2)); 

else 

wgt2=newWgtS€t+newmix[esidx2-newstatebase]; 
float aa.bb; 
float SuinVar[36]; 

mijinmix=(nmixl>ranix2)?nmixl.*njnix2; 
maxiurux=(rmuxl>nmix2)?ruruxl:mnix2; 
if ((ninix!=nmixl) | | (ninix!=ninix2)) ( 
printf ("ex%d and ex%d*s Mixture Number mismatch, they are %d and 
%d\n",esidxl,esidx2,nmixl,nmix2); 

/* Here we only choose the state with max Mixture Number s parameter as the newstate's parameter 

V 

for (int i=0;i<maxnmix;i++) { 
if (nmix2>nmixl) { 
if (esidx2<newstatebase) { 
mean2=hinmparam->Mean(mix2-hi); 
var2=hmnTiparam->Var(mix2+i); 
}else{ 

mean2=mean2_r4'i*DSvecSize; 
var2=var2_r+i*(DSvecSize+l); 

} 

CGLTopl2[i].combWeight=wgt2[il; 
for (int j=0,j <DSvecSize,j++) { 

CGLTopl2[i].combMean[j]=mean2(j]; 

CGLTopl2[il.combVar(jl=var2[j+ll; 

} 

lelse{ 

if (esidxl<newstatebase) { 

meanl=hmmparam->Mean(mixl+i); 

var l=hmmparam-> Var(mixl +i); 
}else{ 

meanl =meanl_r+i*DS vecSize; 
varl=varl_r+i*(DSvecSize+l); 

1 

CGLTopl2[il.combWeight=wgtl(i]; 
for (int j=0,j<DSvecSize,i++) { 

CGLTopl2[i].combMean[j]=meanl[j]; 

CGLTopl2[i}.combVar(jl=varl(j+ll; 

} 

} 

} 

return -1.0; 

1 
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for (int ii=0;ii<nmix;ii+-i-) I 
aa=exp(wgtl[u]); 
al=a->occ*aa; 
if (esidxl<newstatebase) | 
meanl=hmmparam->Mean(mLxl+u); 
var 1 =hmmparam-> Var(mixl +11); 
}else{ 

meanl=meanl_r+ii*DSvecSize; 
varl=varl_r+ii*(DSvecSi2e+l); 

I 

for (int jj=0,jj<nmix,"jj-M-) { 
bb=exp(wgt2[jj]); 
bl=b->occ*bb; 
if (esidx2<newstatebase) { 

meari2=hinmparam->Mean(rrdx2+jj); 

var2=hininparain->Var(rnix2-}-jj); 
}else{ 

mean2=mean2_r+jj*DSvecSize; 
var2=var2_r+jj*(DSvec5ize+l); 

} 

CGL[ii*ninix+jj].coinbCost=0.0; 
if (bhat) { 
for (int i=0;i<DSvecSize;i++) { • 
SumVar[il=:05*(1.0/varl[i+ll+1.0/var2[i+l]); 

} 

} 

for (int i=0;i<DSvecSize;i-M-) { 

CGL[ii*nmix+jj]xombMean[i]=(al*meanl[i]+bl*meari2[i])/(al+bl); 

CGL[ii*nmix+^jj]xombVar[i]=(al^l.O/varl(i+l]+(meanl[i]-CGL[ii*ni^ 
(meanl[il-CGL[ii*ninix+jj]xombMean[iI))+bl*(1.0/var2[i+l]+(mean2[i]<:GL[u^ 
(mean2[i]-CGL[ii*nmix+jj]xombMean[il)))/(al+bl); 

CGL[ii*nmix+jj].combWeight=(al+bl)/(a->occ+b->occ); 

if (!bhat) { 

CGL[ii*nirux+jjlxombCost+=(al+bl)nog(CGL[ii*ninix+jj]xombVar[i])H-ainog(varl[ 
+binog(var2[i+l]); 
) 

/* Bhattacharyya Distance */ 
else { 

CGL[u*nmix+jj]xombCost+=05nos(SumVar[i])+0.25*(los(varl[i4-l])+log(var2[i+l])) 
CGL[ii*rurux+jj].combCost+=0.125*(meanl[il-mean2[i])*(rneanl[i]- 
mean2[il)/SumVar[il-K)5*log(SumVar[i])H<l.25*(log(varl[i+ll)+log(var2[i+l])); 

1 

} 

1 

} 

/* sort the nmix*nmix distances to nmix */ 

for (int toprimix=0;topnmix<nimx;topnmix++) { 

idxl [ topnmix] =nmix+l ; 

idx2[toprunix]=ninix+l; 

1 

bool notfirst=false; 
bool Rep=true; 

for (int toprunix=0;topnrnLx<xm:vbc;topnxriix4H-) { 
float minCost=1.0e+10; 
Rep=true; 

for (int iii=0;iii<ninix;iii+-i-) { 
for (int jjj=0,jjj<nnux,jjj-M-) { 
if (notfirst) { 

for (int counttopn=0;counttopn<topnmix;counttopn+-r) { 
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if (idxl[coimttopn]==iii) Rep=false; 

1 

for (int counttopn=0;counttopn<topnmix;counttopn++) { 
if (idx2(coiinttopnl==jjj) Rep=false; 

} 

if ((CGL[iii*nirdx+jjj].combCost<minCost)&&Rep) { 
rrdnCost=CGL(iii*nnux-i-jjj].combCost; 
idxl[topnmixl=iii;idx2[topnniLx]=jjj; 

} 

Rep=true; 

} 

else { 

if ((CGL[Lii*nmix+jjj]xombCost<minCost)) { 
iTunCost=CGL[iii*miux^-jjj].combCost; 
idxl [ topninix]=iii;idx2[ topnimx]=jjj; 
//notfirst=true; 

} 

I 

} 

1 

notfirst=true; 

CGLTopl2[topninix]xombCost=minCost; 

for (int ttt=0;ttt<DSvecSize;m-M-) { 
CGLTopl2[topninix]xombMean[ttt]=CGL[idxl[topimux]*nmix+idx2[topn^^ 
CGLTopl2[topninix]xoinbVar[ttt]=1.0/CGL[idxl[topnirux]*nnu^ 

CGLTopl2[topmiux]xombWeight=log(CGL[idxl[topmnixl*rm 

} 

. CGLTopl2[nmixI.combCost=0.0; 
for (int i=0; i<ninix; i-M-) 
CGLTopl2[ninix].combCost+=CGLTopl2[i].combCost; 

return CGLTopl2[iunix]xoinbCost; 

} 

void DSizeModel::CopyParaml(aDcsTreeNode *nl, aDcsTreeNode *n2, int newstate) { 
int esidxl=nl->state; 
int esidx2=n2->state; 
int rm^Lxl4imix2,maxninix; 
int *newTunix=newNmixSet; 
if (esidxl<newstatebase) { 
ninixl=hinmparam->ExState(esidxl)->ninix; 

I else { 

nmixl=newnmix[esidxl-newstatebase]; 

} 

if (esidx2<newstatebase) { 
nmix2=hixunparam->ExState(esidx2)->niriix; 

} else { 

ninix2=newninix[esidx2-newstatebase]; 

) 

maxnniix=(nniixl>riiriix2)?iunixl:nmix2; 
/ / printf ("nmaxmix=%d\n",nmaxinix); 
int *newmix; 
if (newstate==0) ( 

newTnix=newMixS€t; 

ne wmix [news ta te] =0; 
}else{ 

newmix=newMixS€t; 

newTnix[newstatel=newmix(newstate-l]+newTunix[newstate-l]; . 

1 

newnmix[newstatel=maxnmix; 
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float *newmeai\=newMeartSet+newmbc[newstateI*E)SvecSize; 
float *newvar=newVarSet+newmix[newstatel*(DSvecSize+l); 
float ♦newwgt=newWgtSet+newmix(newstatel; 
for (int i=0;i<maxnxnix;i-M-) | 
ne\Nrwg;t[i]=FinalCGLTopl2[i].combWeight; 
float tempvar=0.0; 
for (int j=0,j<DSvecSize,i++) { 
newmean(jl=FinalCGLTopl2(i].combMean[j]; 
newvar(j+l]=FinalCGLTopl2[i].combVar(j]; 
tempvar+=log(2.0*M_PI)-{float)log((double)newvar(j-i-l]); 

1 

newvar[01=teinpvar; 

newinean+=DSvecSize; 

newvar+=(DSvecSize+l); 

1 

} 

int DSizeModel::GaussianDistance_voc(void) { 
printfC entering \n"); 

printf( newstatebase=%d\n",newstatebase); 
int newstate=0; 

aDcsTreeNode *temp_node=NULL; 

FinalCGLTopl2=(aCombGaiiss*)mem->MaUoc(sizeof(aCombGauss)*nmix,true); 
CGLTopl2=(aCombGauss ♦)mem->Malloc(sizeof(aCombGauss)*(nmix+l),true); 
TempCGLTopl2=(aCombGauss *)mem->Malloc(sizeof(aCombGaiiss)*(ninix+l),true); 
CGL=(aCombGaxiss*)rnem->MaUoc(sizeof(aCombGauss)frurux*iimix,true); 
idxl=(int *)mem->Malloc(sizeof(int)*ninix,true); . 
. idx2=(int *)meni->Malloc(sizeof(int)*ninix,true); 
for (int i=0;i<nmix+l;i++) { 
if (i<nmix) { 

FinalCGLTopl2[ilxombMean=(float*)mem->Malloc(sizeof(float)*DSvecSi2e,true); 
FinaiCGLTopl2[i]xombVar=(float*)mem->MaIloc(sizeof(float)*DSvecSize,true); 

CGLTopl2[i]xombMean=(float*)mem->Manoc(sizeof(float)*r:)SvecSize,true); 
CGLTopl2[i].combVar=(float *)mem->Malloc(si2eof(float)*DSvecSize,true); 
TempCGLTopl2[il.combMean=(float*)mem->MaUoc(sizeof(float)*DSvecSize,true); 
TempCGLTopl2[i] .combVar=(float *)mem->Malloc(sizeof (float)*DSvecSi2e,true); 

} 

for (int i=0,i<nmix*nmix;i++) { 
CGL[i].combMean=(float *)mem->MaJloc(sizeof(float)*DSvecSize,true); 
CGL[il.coir\bVar=:(float*)inem->Malloc(sizeof(float)*DSvecSize,true); 

} 

for (int i=04<mona->monoFhnNum;i++) { 
if (i==DSte€) continue; 
for (intj=0,j<3,j++) { 
if(DStrees[i][j]){ 

/* normal models*/ jv ■ . ^ x 

printf("Now we are processing %d exstate s %dth state, start newstate=%d\n ,i,j,newstate); 
newstate=VocMergeNodes>I(&DStrees[il[j]->root,++DStrees[i][j]->root.magic,^^ 
temp.node=VocMergeiModes.Refine(&DStrees[il[j]->root,-rH-DStrees[i][j]- 

>rootmag;ic,temp_node); 
} 

} 

} 

mem->Free(CGL); 
mem->Fre€(FinalCGLTop 12); 
mem->Free(CGLTop 12); 
mem->Free(TempCGLTop 12); 
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mein-> Free( idx 1 ); 
mem->Free(idx2); 
return 1; 

I 

DSizeModel::DSizeModel() {}. 

void DSizeModel::Initialize(int tee, DcsTree ***trees, int vecSize, int hmmMix, aHM\rDcsTreeNode 
*dnodes, 

int nodeNumi, bool isbhat) { 
nmix=hinmMix; 
DStee=tee; 
bhat=isbhat; 
DSvecSize=vecSize; 
DStrees= trees; 
DSdcsNodes=dnodes; 
DSnodeNum=nodeNuin; 

1 

Thus, a method and system to scale a decision tree-based hidden markov model 
(HMM) for speech recognition. The method and system described are particularly 
suitable for automatic speech recognition systems (ASR). In the foregoing • ' 
specification, the invention has been described with reference to specific exemplary^ 
embodiments thereof. It will, however, be evident that various modifications and 
changes may be made thereto without departing from the broader spirit and scope of 
the invention as set forth in the appended dainns. The specification and drawings are, 
accordingly, to be regarded in an illustrative sense rather than a restrictive sense. 
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